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I. INTRODUCTION 



There is a growing evidence that Universe is undergo- 
ing accelerated expansion at present. The observations 
related to supernova, CMB and galaxy clustering all to- 
gether strongly point towards the compelling possibil- 
ity of late time accelerated expansion of Universe. The 
acceleration of Universe can be accounted for either by 
modifying the left hand side of Einstein equations or by 
supplementing the energy momentum tensor by an ex- 
otic matter with negative pressure, popularly known as 
dark energy. In past few years there have been tremen- 
dous efforts in modeling the dark energy. They include 
scalar field models, some models of brane worlds and spe- 
cific compactification schemes in string theory which have 
been shown to mimic the dark energy like behaviour. A 
wide variety of scalar field models have been conjectured 
for this purpose including quintessence 0, K-essence0, 
tachyonic scalar fields with the last one being originally 
motivated by string theoretical ideas. All these models 
of scalar field lead to the equation of state parameter w 
greater than or equal to minus one. However, the recent 
observations do not seem to exclude values of this param- 
eter less than minus one'57^. It is therefore important to 
look for theoretical possibilities to describe dark energy 
with w < —\ call ed^ hantom en er^^j^^^^|^[^^| ^ 

Specific models in brane world or non-minimally coupled 
scalar fields may lead to phantom energy [i^, . In our 
opinion, the simplest alternative is provided by a phan- 
tom field, scalar field with negative kinetic energy. Such 
a field can be motivated from S-brane constructs in string 
theory m [H H HI El EHil El . Historically, phan- 
tom fields were first introduced in Hoyle's version of the 
Steady State Theory. In adherence to the Perfect Cos- 
mological Principle, a creation field (C-field) was for the 



first time introduced [50| to reconcile with homogeneous 
density by creation of new matter in the voids caused 
by the expansion of Universe. It was further refined and 
reformulated in the Hoyle and Narlikar theory of gravita- 
tion (see also Ref.[53| on the similar theme). Though 
the quantum theory of phantom fields is problematic |53| . 
it is nevertheless interesting to examine the cosmological 
consequences of these fields at the classical level. 

Models with constant w < —1 lead to unwanted future 
singularity called hig rip^]. This singularity is char- 
acterized by the divergence of scale factor after a finite 
interval of time. Infact there is no profound reason to 
assume that the equation of state parameter w would re- 
main constant in future if it has been changing all the 
way since the end of inflation to date. Keeping this in 
mind, specific scalar field models were proposed to avoid 
the cosmic doomsday [Tsl . It requires a special class of 
phantom field potentials with local maximum. In this pa- 
per we examine the nature of future evolution of Universe 
with a phantom field evolving in a more general class of 
potentials. The future course of evolution is shown to 
critically dependent on the steepness of the underlying 
potential. 



II. COSMOLOGICAL DYNAMICS IN THE 
PRESENCE OF PHANTOM FIELD 

The Lagrangian of the phantom field minimally cou- 
pled to gravity and matter sources is given bv |20ll2q 



(1) 



where Csource is the remaining source term (matter, ra- 
diation) and V{4)) is the phantom potential. The kinetic 
energy term of the phantom field in JQ) enters with the 
opposite sign in contrast to the ordinary scalar field (we 
employ the metric signature, -,-|-,-|-,-f ). In a spatially fiat 
FRW cosmology, the stress tensor that follows from 
acquires the diagonal form — diag {—p,p,p,p) where 
the pressure and energy density of field cf) are given by 
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The corresponding equation of state parameter is now 
given by 



w 



P0 f - ■ 



(3) 



For p0 > 0, w < —1. 

The equations of motion which follow from are 



H = - 
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where the background energy density due to matter and 
radiation is given by 



Pb - — + — 5" 
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As in the previous work |lfi[l2()| we start with the phan- 
tom field energy density being dominated by the density 
of the ordinary matter. In this situation, an initial kinetic 
term of the field rapidly decreases due to the friction term 
2>H(j) in ijHJl and as a result the phantom field "freezes" 
at some point waiting for the moment its energy den- 
sity becomes comparable with the matter (see figure^. 
Subsequently, the field dynamics "switches on" and the 
nature of future evolution then depends on the shape of 
the phantom field potential V[<f)). 

The case of potentials with a local maximum was al- 
ready described in [THl l2fil | . The phantom field asymp- 
totically reaches the maximum with a possible oscillating 
transient stage. The future asymptotic regime is the De- 
Sitter one with w — 1 and A = Vmax/M"^. 

Potentials other than these can lead to a variety of 
future evolutionary regimes. We begin with a simplest 
potential of a massive phantom field 



2 J,2 



V{(j}) = m^' 

When the phantom energy begins to dominate and 
small, we have, see equation ^ and 1© 
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Once such a regime is established its stability is guar- 
anteed (see figure 2). Indeed, the ratio 4>'^/2V of the 
phantom field is inversely proportional to (jP and tends 
to zero making the kinetic term subdominant forever. 
This regime is analogous to the slow-roll regime for a 
normal scalar field and can be named as " slow-climb" . 
The effective equation of state tends to the De-Sitter one 



(see figure 3) with a slow growing energy density (figure 
1). As a result, the final state of such Universe is for- 
mally different from both the De-Sitter and the "big rip" 
outcome - it would take infinite time to reach an infinite 
energy density in this case. However, such Universe will 
steadily reach the Planck density for the finite time and 
the classical physics breaks down here. 

It can easily be checked that the condition (fp' /2V 
is satisfied for any potential, which asymptotically re- 
duces to power-law. Indeed, if the kinetic term is sub- 
dominant, equation ^ gives 



In the regime of phantom field dominance we have 
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It is clear from H12|l that (fp' /2V is proportional to (/)~^ 
for any power-law potential V{4>) ~ 0". 

However, the formally "slow climb" regime for steeper 
power-law potentials becomes the regime of a fast growth 
of the phantom field. The potential V = Xcf)^ leads to 



H : 



corresponding to 



3M2 



V = AXcjP 



(13) 
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which results in exponential growth of (j). A bigger power 
index a leads to a "classical big rip" when a formally 
infinite energy density is reached during a finite time. 
Indeed, the general "slow-climb" equation for the phan- 
tom field with a power-law potential 



.a/2-1 



in case of a > 4 leads to 



(15) 



(16) 



For completeness we note that asymptotically fiat poten- 
tials, hke V{(l)) — A{1~ exp (— c0^)) also lead to a "slow 
climb" but with asymptotically De-Sitter like behavior 
(A ^ F(oo)/M2). 

The approximation (fp /2V breaks down for ex- 
ponential and steeper potentials. In case of exponential 
potentials[23, Hao and Li found an attractor solution 
with w tending to a constant value less then —1 making 
the "big rip" inevitable. 
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FIG. 1: Evolution of energy density is plotted against the 
scale factor. Solid line corresponds to the phantom field en- 
ergy density in case of the model described by the potential: 



with m ~ 10""" M„. The dashed and dotted 



lines correspond to energy density of radiation and matter. 
Initially, the phantom field mimics the cosmological constant 
like behaviour and its energy density is extremely subdom- 
inant to the background (the system is numerically evolved 
starting from the radiation domination era with pr = IMeV^) 
and remains to be so for most of the period of evolution. The 
phantom field (j) continues in the state with w = —1 till the 
moment approaches pj. The background ceases now to 
play the leading role (becomes subdominant) and the phan- 
tom field takes over and starts climbing up the potential fast. 
At late times, the field energy density catches up with the 
background, overtakes it and starts growing {w < —1) and 
drives the current (the value of the scale factor a ~ 4 x 10^ 
corresponds to the present epoch) accelerated expansion of 
the Universe with il^ ~ 0.7 and ilm ~ 0.3. Initial value of 
the field was tuned to get the present values of il. The field 
then enters into the slow climbing regime allowing the slow- 
ing down of growth of and making w —> —1 asymptotically 
which is an attractor in this case. 



When we turn to potentials steeper than an exponen- 
tial, like V{(f)) = exp (ccf)'^) or potentials with an infinite 
potential wall, we get another type of a future singular- 
ity. It is evident from equation H12() that the kinetic term 
grows faster than the potential one, and at a particular 
epoch becomes important. However, p^, — V{(j)) — 0^/2 
always grows as the equation Q in the epoch of phan- 
tom dominance prevents H and correspondingly p^ (see 
equation 0) from decreasing. As a result, both poten- 
tial and kinetic terms increase rapidly with less rapid 
increase of p^. Our numerical studies indicate that the 
ratio of kinetic and potential terms tends to — 1 and the 
state equation parameter u blows up to — oo at some fi- 
nite value of the scale factor a (figure 3). Potential and 
kinetic terms of the phantom field also become infinite at 
this moment. This type of cosmological singularity, when 
equations of motion become singular at a finite value of 
scale factor are already known in some models of brane 
worldsjsil, tachyonic field[3^ and Gauss-Bonnet cosmol- 



FIG. 2: Display of the phase portrait (plot of Y2 = 4>/Mp ver- 
sus Yi = (j)/Mp) for the phantom field in case of the quadratic 
potential V{(l>) = m^(^^. The fi gure shows that the trajec- 
tories starting anywhere in the phase space move towards a 
configuration(s/ow) climb regime) with const asymptoti- 
cally(the convenient choice of mass parameter in the potential 
used here corresponds to the value of the const equal to 2/ \/3) • 
This picture is drawn in absence of the background energy 
density and allows to probe a wider class of initial conditions; 
the asymptotic regime is independent of the background. 
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FIG. 3: Equation of state parameter w for the phantom field 
is plotted against the scale factor. The solid line corresponds 
to the exponential potential V{(j)) ~ e^'i'/'^'p^ the dotted line 

to F(0) ~ e^"*'/*^p and the dashed line to V{(l>) ~ 0^ The ex- 
ponential potential exhibits a critical behaviour in which case 
w fast evolves to a constant value less than minus one (nu- 
merical value depends upon A which was taken to be equal to 
one for convenience) . For potential steeper than the exponen- 
tial w fast evolves towards larger and larger negative values 
(dotted line) leading to a singularity at some finite value of 
the scale factor. In case of a less steeper potential than the 
exponential, the equation of state parameter w first decreases 
fast and then subsequently turns back towards minus one and 
approaches it asymptotically (dashed line). 
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ogy |5g. However, we should note that all these scenarios 
lead to divergence in the second derivatives of metrics. 
The matter content of the Universe and Hubble parame- 
ter remain finite in these models. In case of the phantom 
field, the energy density itself blows up at some finite 
value of a. 



III. CONCLUSIONS 

We have studied the possible future regimes of Uni- 
verse with a phantom field. The opposite sign in the 
kinetic term pushes the field towards maximum of its 
potential. This feature leads to a cosmological singular- 
ity in the future for phantom field potentials unbounded 
from above. The nature of this singularity crucially de- 
pends upon steepness of the potential. The analysis pre- 
sented above allows us to classify the future cosmological 
behavior in the following types: 

• For asymptotically power-law potentials U(0) (p" 
with a < 4 we obtain lo ^ —1 and — > oo for 
t — > oo. 



• For a > 4 though w — s- — 1, but the infinite energy 
density is reached during a finite time, referred to 
as big rip. 

• For exponential potentials, uj tends to a constant 
value less than —1 again leading to a big rip singu- 
larity. 

• Potentials steeper than exponential lead to a sin- 
gularity with If ^ — oo for a finite value of scale 
factor. 
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